Specification 

Be It Known That WE, MICHAEL A. METCALF a citizen of 

the United States of America, resident of San Diego, County 

of San Diego, State of California, and JOHN M, LAND citizen 

of the United Kingdom, resident of Godalming, Surrey, 

5 England have invented a new and useful 

METHOD AND APPARATUS FOR COMBINED MEASUREMENTS OF 
CONCENTRATION, DISTRIBUTION AND FLOW VELOCITY OF SUSPENDED SOLIDS 

of which the following is a specification: 
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Field of the Invention 

This invention relates to measurement of water quality, 
and more specifically to instruments used for measuring the 
concentration of suspended solids in a liquid flowing 
5 through a conduit . 

Background of the Invention 

The safe disposal of sewer effluents and water 
treatment requires measurements of suspended solid 
10 concentration in treatment pipes, culverts and other 

conduits . 

A sewer effluent containing a low value of suspended 
solids may be subject to a reduced treatment or no treatment 
at all. A plant treating sewer effluents from several 
^ municipalities may bill those municipalities in accordance 

with the volume flow and the total suspended solids 
concentration measured on the intake conduits from each 
municipality . 

In the past, intrusive periodical samplings of solids - 
2q carrying water were taken and analyzed in laboratories to 

provide needed information about the solids content of the 
water. 

Various methods have been proposed to electronically 
measure total solids concentration in a volume of liquid, 
25 but no instrumentation has been provided for concurrently 

measuring suspended solids concentration as well as flow on 
a continuous, real time base. 
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This invention results from an attempt to provide a 
more practical instrumentation for use in sewer and water 
collection systems, distribution systems and treatment 
plants . 

5 

Summary of the Invention 

The principal and secondary objects of this invention 
are to provide a method and apparatus for continuously 
measuring the total solids concentration in a liquid flowing 

10 through a pipe or conduit as well as the speed of various 

layers of flow within that pipe and calculating flow rate. 

These and other valuable objects are achieved by 
placing an acoustical transmitter on the inside wall of a 
pipe or other conduit. The transmitter emits two pairs of 

-15 obliquely divergent beams, one beam of each pair is aimed 

downstream from the other one in the pair. Echo signals of 
the emitted waveform are divided into sampling periods 
corresponding to discrete volumes of the liquid in the pipe 
distributed along each beam. Doppler frequency shifts 

2q received from the beams are interpreted into velocity 

measurements of the suspended solids assumed to be flowing 
ambiently with the water within a number of flow layers 
within the pipe. The intensity of the backscattered echo 
signals are translated into solids concentration values 

25 using in sita calibration measurements and simplified 

algorithms, combined with iterative routines that adjust 
calibration parameters using data obtained from previous 
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measurements on an adjacent layer in a continuous self- 
correcting process. 

Brief Description of the Drawing 

Figure 1 is a diagrammatical perspective view of a 
conduit equipped with a flow and concentration measuring 
transducer; 

Figure 2 is a diagrammatical cross- sectional view of 
the conduit of Figure 1 taken along line 2-2 thereof; 

Figure 3 is a diagrammatical cross-sectional view of 
the conduit of Figure 1 taken along line 3-3; 

Figure 4 is a block diagram of the concentration, 
distribution and flow velocity measurement instrument; 

Figure 5 is a block diagram of the electronic 
processing unit; and 

Figure 6 is a block diagram of the transducer 
circuitry. 

Description of the Preferred Embodiment of the Invention 

Referring now to the drawing, there is shown in Figures 
1-3, a suspended solids velocity and concentration 
transducer 1 installed at the bottom of a pipe 2 carrying 
water 3 containing suspended solids 4 and flowing in the 
direction indicated by the arrow 5 in Figure 1. 
Orientations within the pipe can be specified by reference 
to the relative horizontal axis H which would be the central 
major axis of the generally cylindrical pipe and relative 
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vertical axis V which is perpendicular to H and 
substantially normal to the surface of the water. The plane 
P WF is aligned vertically with the flow and contains both V 
and H axes. The plane P CF is aligned vertically across the 

5 flow containing the V axis and perpendicular to the H axis. 

Piezoelectric ceramics in the transducer emit an acoustical 
waveform consisting of short pulses along four narrow beams 
6, 1, 8, 9 pointing in different directions. The beams are 
grouped into two pairs . 

10 The first pair 6/ 7 has a first upstream beam 6 

oriented at a first angle Lg cf from the cross-flow plane P CF 
tending into the flow and a second angle Lg wf from the with- 
flow plane P WF tending toward the side wall of the pipe. 

The first angle Lg c f is preferably between about 5 and 

15 80 degrees, more preferably 10 and 30 degrees, and most 

preferably about 20 degrees. The second angle Lg w f is 
preferably between about 0 and 60- degrees, more preferably 
between about 0 and 30 degrees, and most preferably about 10 
degrees . 

20 The first pair has a second down stream beam 7 oriented 

in a mirror image of the first beam 6 about the cross- flow 
plane P CF - 

The second pair 8, 9 has a first upstream beam 8 
oriented at a first angle Lg cf from the cross-flow plane P CF 
25 tending into the flow and a second angle L Qwf from the with- 

flow plane P WF tending toward the side wall of the pipe. 

The first angle L 8cf is preferably between about 5 and 
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80 degrees, more preferably 20 and 6 0 degrees and most 
preferably about 3 0 degrees. The second angle I^swf is 
preferably between about 0 and 60 degrees, more preferably 
between about 0 and 30 degrees and most preferably about 10 
degrees . 

The second pair has a second down stream beam 9 
oriented in a mirror image of the first beam 8 about the 
cross flow plane P^F" 

As more specifically shown in Figure 2, each pair of 
beams preferably lie within a plane which is aimed at a 
transversal angle from the vertical and from the other pair 
of beams. Echo signals of the pulses are backscattered from 
the suspended solids 4. Since these solids have motion 
relative to the transducer, the echo signals are Doppler- 
shifted in frequency. A fifth ceramic transducer 10 mounted 
in the center of the transducer and aimed vertically is used 
to measure the depth of flow. 

As illustrated in Figures 4 and 5, the system includes 
an electronic processing unit 11 that receives an input 
signal RCVS-SIG from the transducer 1 and converts it into 
digital data D-OUT that is fed to a data processor 12. The 
system divides the echo signals into discrete regular 
intervals for samplings that correspond to different, 
discrete volumes at different ranges of the flowing liquid. 
Velocity is calculated from the frequency shift measured in 
each sampling. The result is a profile of linear 

distribution of velocity along the beams illustrated by 



velocity profile lines 13 in Figure 2. Each of the small 
divisions 14 shown along beam 9 of Figure 2 represent an 
individual velocity measurement in a discrete volume known 
as a depth cell or bin. The velocity profile lines 13 are 
generated from velocity data measured by the upstream and 
downstream beam of each pair. 

Since Doppler measurements are directional, only the 
component velocity along the direction of transmit and 
receive is measured. Narrow acoustic beams are used to 
accurately determine the horizontal velocity of the flow. 
The accuracy of the measurement is also enhanced by using 
gating times for the samplings that correspond to small 
volumes of a minimum size of approximately 5 centimeters in 
length and 4 centimeters in diameter and approximately 1/25 
of the beam length to conduit wall. Potential bias in the 
return energy spectrum due to range dependent variables is 
avoided. The result is a very precise measurement of the 
vertical and transversal distribution of flow velocities. 
The velocity data from the two pairs of beams are entered 
into an algorithm to determine a mathematical discription of 
the flow velocities throughout the entire cross-section of 
the liquid. The algorithm fits the basic funtions of a 
parametric model to the actual data. The results predicts 
flow velocity at all points throughout the liquid. These 
results are integrated over the cross-sectional area to 
determine the discharge. 

The key benefit of this approach is that the system 



will operate accurately under different hydraulic 
conditions. As hydraulic conditions change, the change will 
manifest itself in the distribution of velocity throughout 
the depth of flow. As the system is measuring the velocity 
distribution directly, it will adapt to the changes in 
hydraulics, and generate a flow pattern that is 
representative of the new hydraulic conditions, insuring an 
accurate calculation of the flow rate. 

The measurement of distribution and concentration of 
the suspended solids 4 stems from a requirement of the 
frequency measurement circuits within the electronic 
processing unit 11 to receive a near-constant input voltage. 
An acoustic pulse emitted from the transducer 1 has a 
certain initial intensity that progressively diminishes as 
it travels through the liquid and is scattered from 
suspended particles. The reflected energy detected by the 
transducer is a very small fraction of that which was 
emitted. In addition, the return intensity of the 

backscattered echoes vary considerably according to, but 
not limited to, the range to the point of energy reflection, 
the concentration of reflecting particles, and the water 
temperature. The received weakened signal RCV-SIG passes 
through an amplifier RCV-AMP that brings it to the level 
required by the frequency measurement circuits. A large 
loss of signal strength requires a large degree of 
amplification. The degree of amplification which is 
required is thus a measure of the loss of signal strength 
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and, inversely, a measure of the intensity of the 
backscatter. The amount of amplification required is 
provided by the Receive Signal Strength Indicator RSSI. It 
is this measurement which facilitates the estimation of 

5 suspended solid concentration in the water column. In other 

words, the intensity value of the backscattered signal is 
translated into concentration values of the suspended solids 
using the data processor 12 and the frequency shifts of the 
receive signal are interpreted as an indication of velocity 

10 of the solids in the flow of liquid. 

The transducer output signal RCV-SIG is routed to a 
receive coupling transformer RCV-XFMR that provides 
isolation and impedance matching. The signal is further 
amplified and bandwidth limited by the high gain selective 

15 log amplifier RCV-AMP. The amplified receive signal is fed 

to a pair of frequency mixers FMIX, where the signal is 
mixed replicatively with the local oscillator signal 
frequency LO. The desired base band signal, which is the 
difference frequency of the receive signal and the local 

2o oscillator frequency, is obtained by passing the mixer 

output signals through a pair of low pass filters LPF. The 
base band signals contain now the entire Doppler spectrum 
without the carrier signal. The mixers, are quadrature 
mixers, where an in-phase and a quadrature signal are 

2 5 obtained. Both signals are needed for the correlator CORR, 

which performs the basic digital signal processing. The 
in-phase and quadrature signals are buffered by a pair of 
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first-in/first-out buffers FIFO and become part of the data 
D-OUT output to the data processor 12. The echo signals are 
also fed to a low-pass filter RSSI-LPF, then to a buffer 
RSSI-BUFF and digitized by an analog- to-digital converter 
RSSI-ADC. Finally, they are fed to the data processor as 
part of the output signal D-OUT. 

A timing generator TMG-GEN generates all signals needed 
for the transmitter and receiver, such as the transmit 
signals, transmit enable, and the local oscillator 
quadrature signal for the mixer. The frequency of the 
pulses emitted by the transducer is selected to maximize 
resolution while maintaining adequate profiling range. In 
the present, a frequency of about 1.23 MHz is preferred. A 
transmitter amplifier XMT-AMP acts as a power driver which 
buffers the logic level signal generated by the timing 
generator, and drives a transmitter output transformer XMT- 
XFMR, the transmit transformer also provides isolation 
between the electronic processing unit and the transducer. 

The transmit current is monitored by a current 
transformer CURR-XFMR. Its output signal is scaled and 
digitized by an analog- to-digital converter CURR-ADC and is 
used as part of a built-in self-test by the data processor. 
All timing generator setups are fully programmable, and are 
down-loaded by the data processor to the timing generator's 
on an internal RAM. 

The data processor is able to read back, the timing 
setup data and the digitized current sensed data. 
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The transducer illustrated in Figure 6 includes a 
signal transmit and receive multiplexer controlled by the 
electronic processing unit. The transmit signal generated 
by the timing generator TMG-GEN of the electronic processing 
unit passes through a coupling transformer XFMR which 
provides isolation and impedance matching. The signal is 
then multiplexed to one of the five beams by a selectable 
transmit/receive switch T/R-SW6-T/R-SW10 . At the end of the 
transmitting phase, the multiplexer is deselected. As the 
echo signal is received from the Piezoelectric ceramics P-6 
through P-10, it passes through the switch and one of the 
selected band pass filters BPF-6 through BPF-10 then to a 
preamplifier PA- 6 through PA- 10 and is converted to a 
differential signal by a wide band signal transformer XFMR 
15 before being fed to the electronic control unit. A 

thermiistor T- SENSE is used to measure the transducer's 
ambient temperature. The temperature signal is scaled and 
buffered by an amplifier S-AMP, digitized by the transducer 
control XDCR-CTL and fed to the electronic control unit. 
20 The following simplified formula of the acoustic theory 

governing backscatter from particles suspended in the water 
column identifies the main factors that contribute to the 
determination of suspended solids concentration from the 
intensity data. 



25 



E = SL + SV + Constant - 20\o°(R) - 2a w R 
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where: E 



echo intensity, 
SL s transmitted power, 

SV = backscatter intensity due to the particles 
suspended in the water column, 

a coefficient describing the absorption of energy bv 
the water, s/ y 

R s the distance from the transducer to the 
measurement bin. 



a w 



s a 



The echo intensity E measured by the system i 
relative intensity, coming f rom direct measurement of ^ 

pressure amplitude of the return signal. Although the 
system will clearly recognize variations Qf ^^^^^ 

it cannot determine the exact - 

exact amount of backscatter 

intensity that is due solelv t-n f ho 

soxeiy to the presence of solids. 

Other factors that contribute t-o fh Q , . 

oncrioute to the final intensity must be 

removed . 

The term 201og R i s . simple geometrlc functlon ^ 
account for the spherical spreading of the beam. This term 
can be further refined to include a near-field correction to 
the assumption of spherical spreading. Accurate knowledge 
of the velocity of sound is essential for determination of 
the range to a given measurement bin. The system computes a 
range based on the measured time and a sound velocity. Sound 
velocity is computed by using a user-defined salinity 
(assumed to be o for the sewer environment, and the 
temperature measured at the sensor head. This speed of 
sound is assumed by the method described here.to be constant 
throughout the water column. This is a safe assumption in ' 
sewer environments as the depth of the flow is Umited and 
the flow is well mixed. 

It is not possible to make direct comparisons between 



the measurements made by different instruments unless the 
instruments have been calibrated in the field or the 
laboratory to establish their performance characteristics. 
The terms SL and Constant account for these differences and 

5 also deal with the characteristics of the suspended load. 

Understanding the manner in which instrument performance 
characteristics affect the data and correcting for these 
differences through field calibrations is key to the 
measurement method. This method makes use of the measured 

10 relative backscattered intensities and uses site-specific 

calibrations to compensate for such variations and determine 
a concentration of suspended solids. 

The final two terms ^ and SV, refer to the absorption 
of acoustic energy by the water and the relative backscatter 

15 intensity. Acoustic energy is absorbed by water as it 

passes through it, and ^ is a measure of the amount of 
energy lost in this process. SV is the term we are 
interested in - the amount of energy that is backscattered 
by the presence of solids in the water column. An increase 

20 or decrease in the amount of solids will affect the value of 

SV. An additional term not included in Equation 1 must also 
be accounted for. This term describes the attenuation of 
the acoustic signal due to scattering and absorption by the 
suspended load. 

25 To actually derive the mass concentration of total 

suspended solids TSS in the water from the measured 
intensity, one must determine the true backscattered sound 
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intensity due to the presence of the solids. It can be 
shown that the relationship between the two parameters is as 



follows : 



M(r) = (K< Pma > rf lEilPi e 4r (aw + asJ 

<f> 2 



(2) 



V " nere M(r) " mas s concentration per unit volume at 

range r 

K ~ a constant for the acoustic system 
Prms - backscattered pressure 
a s = particle radius 

<f> = Particle density 



Particle form function 
water attenuation coefficient 



To 
a 



10 s so,,ds attenuation coefficient 

In this equation, the expected mass concentration, M(r) 
is a function of the sediment attenuation coefficient a s , 
which defines how the presence of solids attenuates the 
return signal. Both of these parameters are unknown. 
15 provide a solution to this, it is necessary to use 

numerical method. Computation of the mass concentration is 
performed without using a value of a,. The resulting value 
of M(r) is used to compute a value for ... This computation 
process is iterated in the bin nearest to the sensor to 
derive final values of both sediment attenuation and mass 
concentration. This method is then applied in a step-wise 
manner along the successive bins of the profile. 

The method proposed here to allow for "the measurement 
of solids concentration on site, accomplishes this by "using' 
a simplified version of the above-expression 
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dB is now the measured relative backscatter intensity, 
corrected for spherical spreading and any applicable near 
field effects. S is the relative backscatter coefficient 
that defines the relationship between solids concentration 
and particle size. K s is the site and sensor specific 
constant that corrects for the individual characteristics of 
a specific instrument at a given site. The other terms are 
as in Equation 2 and and represent intensity attenuation by 
the presence of solids (a s ) and water absorption (a w ) . 

The amount of sound energy which is lost, or 
attenuated, due to the absorption by water has been found to 
depend on the frequency of the sound waves (i.e., of the 
instrument) , and the salinity and the temperature of the 
water as follows: 



/ 



a * 91500 



where: 



f 
S 



\MSf T f 2.86/ 



(4) 



the water absorption coefficient in Nepers/m 
the instrument frequency, in MHz 
the salinity in ppt. 



The term, ft, is called the Relaxation Frequency and is 
given by the following expression: 



f<-— 1 

f T = 21.9x10^ 273+T) ' (5) 



where: T is the water temperature in degrees Celsius. 

Using these formulae, the attenuation of the signal due 
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to water absorption, per meter, may be computed and hence, 
using the slant length of the bin, the total attenuation of 
sound through each measurement interval, or "bin", can be 
derived. 

When sound passing through the water column strikes 
particles of suspended solids, the energy is attenuated 
both by scattering and absorption by the solid media. The 
degree of scattering depends on the relationship between the 
frequency and the size of the particle. Frequency may be 
expressed in terms of the Wave Number, k, where: 

k = l£ (6) 

where: f = the instrument frequency in Hertz 

V s = the velocity of sound in m/s. 



When the value of the product ka s (where a s is the 
particle radius) is well below 0.5, Rayleigh scattering 
occurs. This is the range normally found during suspended 
sediment measuring conditions. The sediment constant 51, is 
20 given by: 

Ps 



where: Ps 



= the density of the sediment in kg/m 
i£ = a term relating to compressibility and 
density (normal value 0.18) 



The actual sediment attenuation per meter due to 
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scattering, a x (Nepers/m) , by the solids in a region at 
range, r, and with a mass concentration of M r (kg/m 3 ) is 
then given by: 



a, = M r x £*, 



(8) 



Acoustic energy is also absorbed by sediment in the 
water column. The following expression (Urick, 1948) is 
used to determine the sediment attenuation per meter due to 
this absorption, a 2 (Nepers/m) : 



a 2 = 



M r k(c-lf 


s 


2p s 


s 2 + (a + 6) 2 



(9) 



where: 



Pw 



s = 



4P4 



1+ 



kV, 



2y„ 



and: 



Ps 
Pw 
Vw 



the density of the solid particles 

the density of the water 

the kinematic viscosity of the water. 



The above-expressions are valid for values of the 
product ka s , being much less than unity. This is the region 
where attenuation due to scattering becomes negligible and 
viscous absorption begins to dominate. 

The method combines the attenuation due to scattering 
and absorption into a single Sediment Attenuation 
Co-efficient SAC. The SAC can either be defined by the user 
or can be computed by the processing softare using a 
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nominal, or "effective", particle size in combination with 
assumed values for particle specific gravity (2.7) and 
compressibility (0.18). In both cases, the input values are 
inevitably estimates and must be refined by iteration within 
the software's Calibration Module. 

It is almost impossible, in practice, to calculate from 
first principles the sediment attenuation coefficient for a 
natural population of sediment. The approach used in the 
software is a pragmatic approach in which workable values of 
the coefficient are established using actual observational 
data. It therefore, becomes somewhat academic as to whether 
or not the input values of, for example, particle size and 
compressibility are correct. 

The computer program CP controlling the data processor 
carries out all of the computations required to derive 
sediment concentrations from measured backscatter intensity, 
according to the simplified Equation 3 and a value for 
sediment attenuation derived by calibration. The program 
imports the backscatter data (expressed in instrument 
counts) and ancillary information (e.g., temperature, 
salinity, etc.) from raw system data files. The user may 
either input specific site calibrations values (e.g., values 
for S and Ks) on the keyboard or have the program determine 
these values from entered data. The user also input values 
for TSS that were measured by taking samples of water and 
having them analyzed by a laboratory for TSS values . These 
sample are taken concurrently with the system data so that a 
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direct calibration between system-measured values and actual 
values can be determined. 

The program then takes the raw measured data from the 
first valid measurement bin of each transducer. It 
calculates the various required parameters and solves 
Equation 3, using the iterative routine, to find the 
suspended sediment concentration in the first whole bin and 
in the interval between the transducers and the first bin 
(where the system cannot obtain data) . The computed solids 
concentrations and sediment attenuation are then fed into 
the system again in order to make the relevant corrections 
in the next bin, which is then solved using a similar 
iterative process. The procedure is repeated until the last 
whole bin has been solved. The software then moves on to 
the next ensemble of data. 

Following this, the program displays the calculated 
concentrations and compares them to the actual measured 
concentrations. Errors between the two data sets, 

calculated then actual, are also displayed as functions of 
depth and concentration. The user then adjusts the 
different parameters (S, Ks, and the SAC) to improve the 
correlation between data sets. As this is an iterative 
process, and some parameters are initially given seed 
values, the user adjusts the parameters to increase the 
correlation between the two data sets and reduce the errors 
to zero, while observing the results of their actions in 
real-time . 
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The method can be utilized in the sanitary and combined sewer 
environment, and can provide both historical and real-time 
measurement of TSS in addition to flow measurement. This 
will allow users to measure total mass transport of solids 
within sanitary sewer systems as well as the volume 
transport of water. The method described here uses data 
from multiple beams to accomplish this. The method has the 
ability to look at the spatial distribution of solids 
concentration in the depth of flow as well as total solids 
concentration . 



WHAT IS CLAIMED IS: 



